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We consider a chiral baryon-meson model for nucleons and their parity partners in mirror assign¬ 
ment interacting with pions, sigma and omega mesons to describe the liquid-gas transition of nuclear 
matter together with chiral symmetry restoration in the high density phase. Within the mean-held 
approximation the model is known to provide a phenomenologically successful description of the 
nuclear-matter transition. Here, we go beyond this approximation and include mesonic fluctuations 
by means of the functional renormalization group. While these fluctuations do not lead to major 
qualitative changes in the phase diagram of the model, beyond mean-held, one is no-longer free to 
adjust the parameters so as to reproduce the binding energy per nucleon, the nuclear saturation den¬ 
sity, and the nucleon sigma term all at the same time. However, the prediction of a clear hrst-order 
chiral transition at low temperatures inside the high baryon-density phase appears to be robust. 
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I. INTRODUCTION 

Relativistic heavy-ion collisions probe the state of 
strong-interaction matter at finite temperature and finite 
baryon density mm ■ a particularly interesting region in 
the phase diagram is that of large baryon densities. In 
this work we concentrate on cold nuclear matter in the 
vicinity of the liquid-gas transition. Unfortunately the 
sign problem in Lattice QCD prevents a straightforward 
investigation of this region of the phase diagram through 
Monte-Carlo simulations. Methods to circumvent this 
problem are being actively developed. For example, com¬ 
plex Langevin dynamics has recently been successfully 
used for full QCD at finite density albeit on still rather 
small lattices [5]. Strong-coupling techniques can either 
be used together with hopping parameter expansions for 
very heavy quarks or to derive a graph represen¬ 

tation valid also for light quarks but without continuum 
limit E3E|. More direct evidence of a first-order transi¬ 
tion analogous to that of nuclear matter has so far only 
been seen in G 2 -QCD at finite density, a QCD-like theory 
with fermionic baryons but without sign problem which 
can therefore be simulated with standard lattice tech¬ 
niques M- 

Meanwhile, effective models such as the Walecka model 
[T0HT2] have been used since the 1970s and are still be¬ 
ing used to describe the properties of nuclear matter in 
the low temperature region of the QCD phase diagram 
around the liquid-gas transition. In its simplest form 
the Walecka model consists of a single nucleon field cou¬ 
pled to a neutral scalar field and a neutral vector meson 
field. In the mean-field approximation this model shows 
nuclear-matter saturation, when adjusting the model pa¬ 
rameters to realize the nuclear-matter binding energy 
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Eb — 16MeV per nucleon together with a saturation 
density of no = 0.16 fm -3 . On the other hand neutron 
matter remains unbound in the Walecka model. In the 
chiral Walecka model the fundamental degrees of freedom 
are nucleons coupled to the scalar-pseudoscalar sector in 
a chirally invariant way as well as to a neutral vector 
meson. Early on, however, it was noticed that on the 
mean-field level the chiral Walecka model leads to mass¬ 
less Lee-Wick matter, with massless baryonic degrees of 
freedom in the chirally restored phase. Nevertheless, us¬ 
ing a phenomenological parametrization of the effective 
action at zero-temperature and density as the input, it 
has recently received renewed interest. This is because 
it then proved to be very useful for the investigation of 
nuclear and neutron matter at small temperatures both 
in the mean-field approximation m and beyond mm- 

The difficulties with the chiral Walecka model at zero 
temperature are related to the way the nucleon mass 
is generated in this model, predominantly by dynami¬ 
cal chiral symmetry breaking as in quark-meson-models. 
Unlike constituent quarks, however, the nucleons are not 
expected to become massless when chiral symmetry gets 
restored. This is in fact supported by lattice studies 
in which the effects of chiral symmetry restoration have 
been studied by removing low-lying Dirac modes in va¬ 
lence quark propagators CMS]. An alternative model to 
reflect the fact that the nucleon mass stays finite when 
chiral symmetry gets restored as expected with increas¬ 
ing baryon chemical potential beyond the nuclear-matter 
transition must therefore include the nucleons’ parity 
partners Emm, conventionally identified either with the 
IV (1535) or the heavier IV(1650)P] Such models are called 
parity-doublet or mirror models p3H25) . 

Parity-doublet models have provided very promising 


1 Here we use the ;Y( 1535) but for the purpose of our present study 
one might as well use the 7V(1640) with very minor changes. For 
a recent discussion of both assignments see, e.g., m 
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results already at tree-level and in mean-field (MF) ap¬ 
proaches |26ti37j . This was the primary motivation to ex¬ 
amine the effects of including mesonic fluctuations in such 
a model which is done here using the functional renormal¬ 
ization group (FRG). The paper is organized as follows: 
In Sec. [TT] we review the mirror assignment, the parity- 
doublet model and the role of bosonic and fermionic fluc¬ 
tuations in this model. In Sec. lIIIl we describe our results 
in the different approximation schemes at first vanishing 
temperature, and finally also at finite temperature before 
we summarize and conclude in Sec. UY1 


II. PARITY-DOUBLET MODEL AND THE 
NUCLEAR MATTER TRANSITION 

A. Mirror assignment and chirally invariant mass 

First recall the possible representations for baryons un¬ 
der chiral SU(2)l x SU (2 ) R transformations by consider¬ 
ing the product of three (|, 0) ® (0, |) quark representa¬ 
tions where one assigns, without loss, say the first label in 
the product representations (m, n) to an irreducible rep¬ 
resentation of SU{2)l and the second to one of SU(2) R , 
yielding p2Fj . 

((i 0) © (0, §)) © ((|,0) © (0,1)) © ((i, 0) © (0, i))) 

= 5(§,0)©(0,§))©3((1,§)©(§,1)) (1) 

© ((|,0)© (0,|)) . 

For the iso-doublet of proton and neutron one thus of 
course uses the same (|,0) © (0, |) representation as 
for the quarks, for left and right-handed Dirac fermions 
ipi j r = 1: ^ 75 ip, so that the nucleon fields transform just 
like the pair of up and down quark fields. 

If one now considers two baryon species ipi and if> 2 , 
both in a (1,0) © (0, |) representation, there are two 
possibilities: Assigning the first label of the first species 
say to ipij as before, the first label of the second species 
can either refer to ^ 2 ,i or to il ) 2 ,r- The first assignment, 
where both fermions transform exactly in the same way, 
is called the naive assignment whereas the second assign¬ 
ment, where the transformation properties of the second 
left/right-handed fermion species under chiral transfor¬ 
mations are interchanged relative to the first, is referred 
as the mirror assignment. These two assignments are the 
only possibilities in a two-baryon system, provided both 
species belong to a (|, 0) © (0, |) representation. Explic¬ 
itly, they correspond to the transformation properties 

1pl,l -A 1pl,r “A SlRlpl,r 

i>2,1 -A ^Llp2,l, i>2,r -A Cl R 1p 2 , r 

for the naive assignment, and 

1pl,l -A HlV’I.L V’l.r -A O.Rtpl,r 
1p 2 ,r -A Cl L 1p2,r, 4>2,l -A ^Rlp2,l- 


for the mirror assignment, where £Ir,l £ SU(2) R , R - 
Obviously, the generic kinetic term 

£ k in = ( 4 ) 

i— 1,2 

is invariant under chiral transformations, irrespective of 
the assignment. Conversely, for both assignments a con¬ 
ventional Dirac-mass term 

C m = -mtpiipi ( 5 ) 

explicitly breaks chiral symmetry in either case. Hence, it 
is impossible to write down a chirally invariant mass term 
for a single Dirac-fermion species. Interestingly, however, 
this is no longer true for two fermion species in the mirror 
assignment. In this case a mass term of the form 

C m , mirror =m 0 (^ 2 4’l + Vh^) 

. . . . ( 6 ) 

=m 0 (lpl r 4’ll + 4>li4>2r + 4>l r 4>2l + lr) 

remains invariant under the full chiral SU(Nf) R x 
SU(N f ) R . In the naive assignment such a term would 
of course also break chiral symmetry. The particularly 
appealing feature of the mirrror model thus is the possi¬ 
bility of having a chirally invariant local mass term. 


B. Construction of the parity-doublet model 

The parity-doublet (or mirror) model consists of two 
species of mirror-assigned baryons with a chirally invari¬ 
ant mass term of the form of Eq. (Jg]) . Compared to the 
previous subsection they are introduced as Ni = i)>i and 
N 2 = 75^2 so that N 2 has the opposite parity of N\ and 
the eigenvalues of the mass matrix in the chiral limit are 
both mo, rather than ±?uo for that in Eq. ©, see |2Sj. 
Both baryon species are coupled to the scalar/pseudo- 
scalar meson sector (f> = (cr, 7 ?) in an S'0(4)-invariant 
way. The corresponding Euclidean Lagrangian (includ¬ 
ing baryon chemical potential /is and a vector coupling 
to the co-meson) then reads 

C =N t ($ - /is7o + h x (a + i^fn) + i/i„7 M w M ) N x 
+ N 2 {$ — /7b7o + h 2 (cr - i7 5 f7f) + i/i„ 7 M co A( ) N 2 
+ m 0 (N 1 'y 5 N 2 - N2J5N1) + C mes ■ 

(7) 

Note that the mirror baryon N 2 transforms as N 2 —► 
e -i6 7 t jy 2 UIlc j er ax i a i transformations r whereas the 
original baryon transforms as N x -A e lS 7 T N x , which 
requires opposite signs in their Yukawa couplings to the 
pion. Furthermore, we have assumed the same strength 
i h v for the Yukawa coupling of both to the iso-scalar vec¬ 
tor meson which is imaginary here, in order to describe 
a repulsive four-fermion vector interaction. The mesonic 
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part is given by 


where 


Anes = \djd^ + \f^F^ + U^ 2 ^ 2 ) - ca, (8) 

where F^ v = d^ui" — and the mesonic potential at 

tree-level, i.e. in the microscopic bare action at the ultra¬ 
violet cutoff scale A, or in the mean-field approximation 
is of the form 

f/(0W) = - / ^ 2 + ^ 4 + A 6 0 6 + ”|^ 2 , (9) 

with 4> 2 = a 2 + if 2 and parameters p 2 ,X,Xe. A non¬ 
vanishing pion mass is taken into account by means of 
an explicit linear breaking term ccx, i.e. 

V m = U{(I) 2 ,cu 2 )-c<j. (10) 


ocp — 2m 0 -f- 2p~ T A| (j T h 2 (j“, 

P p = -j- - ^{hi - h 2 ) 2 a 2 (Aml + (/ii + h 2 ) 2 (J 2 ). 

Its zeros m± = ipo (at vanishing 3-momentum p = 0 and 
chemical potential p b = 0) define the mean-field baryon 
masses, which are given by 

m± = ^ (±(h! - h 2 )a + \JAmi + ( ft i + (16) 

In particular, in the chirally symmetric case for a = 0 the 
baryons become degenerate with mass Too again. Evalu¬ 
ating the mean-field grand potential for a sharp momen¬ 
tum cutoff leads to 


The vector meson could in principle be consistently in¬ 
cluded as a fluctuating field, for example in the Stueck- 
elberg formalism [5S. 22], but with m u = 782 MeV its 
mass is comparatively large so that it will not have much 
of an impact as a fluctuating meson field in loops. In¬ 
stead, we therefore simply evaluate its expectation value 
on a given solution for the effective potential in a mean- 
field treatment in the following. Due to the repulsive 
nature of the vector interaction, this expectation value 
will be purely imaginary corresponding to a stationary 
phase approximation for complex saddle points. 

Rewriting the Lagrangian in momentum space we find 


/ i 3 

—-^THogS^ + Vkj (17) 

vo ^ 7r) 

with 

Cl q = -2 Nf ^2 I TTTzyy l°g ( cos h 


(27r) 3 


+ Ms 


= - 2 N fJ2 


+ T log 1 + e 


3„ rP± 


d > 
(27Tp 


Ms I 


2 T 


\4 - MsI 


(18) 


I'jf +^B 


+ T log 1 + e 




C = + C 


( 11 ) 


where we have introduced single-quasiparticle energies 


where fk = (^ ) and 

So” 1 = (12) 

— ijS+fr.i(o-+i7 5 7rT)+ifr,„i/i—7157° m 0 75 \ 

— m 0 75 — if+h-z(a—i^ 5 Sf)+ih v ^— P b 1 ° ) 

C. Fermionic and bosonic fluctuations in the 
parity-doublet model 

In the mean-field approximation the inverse baryon 
propagator is given by 

0-1 _ / -if+h 1 cr-p. B m 0 75 \ 

‘-’O — —mo75 -i f’+h 2 a-p,B J 1 1 

where jls = Ms — ih v uio is the baryon chemical poten¬ 
tial shifted by a potentially non-zero value of the zero- 
component of the w-meson field, with A/i^ = ih v uio > 0, 
in the rest-frame of an isotropic thermal medium in which 
the expectation values of the spatial components van¬ 
ish. The determinant required for the computation of 
the grand potential is given by 

det 70-So" 1 = ({po - iMs) 4 + a P (Po ~ i£s f + 0 P ) 2 ( 14 ) 



For mo = 0 one simply has = \Jp 2 + ■ The expres¬ 

sion for the fermionic part of the grand potential then re¬ 
duces to a sum of two usual ones, as e.g., in quark-meson 
models, with masses m + = h\a and ?n_ = h 2 a, respec¬ 
tively. The number of flavors here is Nf = 2 representing 
two iso-doublets, one for the nucleons (p, n) and one for 
their parity partners. 

The mean-field calculation can be improved by in¬ 
cluding mesonic fluctuations, which is done here within 
the framework of the Functional Renormalization Group 
(FRG), a powerful non-perturbative tool in quantum 
field theory and statistical physics, see |40fPH] for QCD- 
related reviews. The central object in the approach pi¬ 
oneered by Wetterich is the effective average action Tk 
which generalizes the effective action F by introducing a 
coarse graining scale up to which quantum fluctuations 
are included. On a technical level this is achieved by 
means of a regulator function Rj. which acts like a mass 
term in the IR. The RG-scale k is taken from the ul¬ 
traviolet (UV) cutoff scale A down to zero and in turn 
the effective average action interpolates between a mi- 
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croscopic bare action at A and the effective action T at 
k = 0. The evolution of the effective average action with 
k is described by an exact 1-loop equation, however in¬ 
volving full field- and scale-dependent propagators, which 
takes the form 


d k T k = | Tr 


,(2) 


+ R k 


-d k R k 


( 20 ) 


( 2 ) 

where T k Ms the second functional derivative of the effec¬ 
tive average action with respect to the fields. The trace 
includes a momentum integration as well as traces over 
all inner indices. 


where one introduces mesonic single-quasiparticle ener¬ 
gies via 


e% = yJk*+2U' k + 4UZ<P, 

e fc = \Jk 2 + 2U' k , 


(26) 


with the notations U k = -£^U k and U k = qSm * 


U k . 


D. Parameter fixing and phenomenology of the 
nuclear matter transition 


To evaluate ( |20| ) one has to specify an Ansatz for the 
effective action. Here we work in the local potential 
approximation (LPA) where the only scale-dependence 
stems from the effective potential, and correspondingly 
we use 


T k = J d 4 x £\u->u k (21) 

with the Lagrangian from Eq. 0 starting with a bare 
potential Ha at fc = A of the form as given in Eq. ^ . The 
flow for the effective potential decomposes into fermionic 
and bosonic contributions, 

d k U k = d k U kt F + d k U ki B ■ (22) 

Using the 3d-analogue of the LPA-optimized regulator 
from Ref. [45], 

Rk,Fip) = -ip- 7 “ l) 6 ^ 2 -P*) ’ ( 23 ) 


the fermionic contributions to the flow of the effective 
potential can be obtained analogous to that of the quark- 
meson-diquark model for two-color QCD in [4f>|, giving 




2 (k 2 + TOq - e^ 2 ) + {h 2 + hl)a 2 




x ^tanh ( j k 2 T B 'j + tanh ^ 


£ k i* B 

2 T 


(24) 


where we have again used the single-quasiparticle ener¬ 
gies e k defined in (19), here at p 2 = k 2 . As in the mean- 
field approximation, for mo = 0 the two fermion species 
decouple and (24 ) reduces to the sum of two quark-meson 
model contributions to the flow of the effective potential. 


The bosonic contribution to the flow of the effective 
potential is identical to the expression in quark-meson 
models [33! and reads, for the 3c?-analogue of the LPA- 
optimized regulator, 


d k U k ,B 


k 4 

12t r 2 




+ — coth 



(25) 


The most relevant degrees of freedom to describe the 
phase diagram of strong-interaction matter near the nu¬ 
clear matter transition are collective mesonic and bary- 
onic excitations. Since there are no baryons in the vac¬ 
uum one should in principal fix the model parameters 
at T = 0 with a value of the baryon chemical potential 
close to the onset of nuclear matter where one has small 
nucleonic excitation energies. Due to the Silver-Blaze 
property one can in principle equally well fix the model 
parameters at T = /is = 0, however. In the LPA this 
is known to introduce a slight artificial /x^-dependence, 
and hence a Silver-Blaze problem, mainly because the 
curvature masses extracted from the mesonic effective 
potential are not exactly the physical ones. Calculating 
mesonic two-point functions in random-phase approxima¬ 
tion or from their own flow equations |46fl50| one observes 
that the physical pole masses can differ quite significantly 
from the curvature masses. In the quark-meson-diquark 
model for two-color QCD or the quark-meson model for 
QCD at finite isospin density one verifies at mean-field 
level that the RPA-pole masses agree with the onset of 
Bose-Einstein condensation of diquarks or charged pi- 
ons, respectively, as they must. In contrast, one then 
deduces that especially the pion curvature mass can de¬ 
viate from this by up to 30% mm- Beyond mean-field, 
pole masses in present truncations for two-point functions 
are typically considerably closer to such onsets than cur¬ 
vature masses as well. To completely resolve this dis¬ 
crepancy, and to reduce the unnaturally large differences 
between curvature and pole masses, one has to include 
wave-function renormalization beyond the LPA Ell¬ 
in the parity-doublet model these differences are much 
smaller, however, because the baryons are much heav¬ 
ier than the quarks in quark-meson models. Even in 
the chirally symmetric regime the mirror baryons have 
their explicit mass too which is of the order of 800 MeV 
and implies that their fluctuations are suppressed below 
RG scales k ~ 800 MeV. A manifestation of this is that 
we have to start the flow already in the broken phase 
when using a typical UV cutoff scale A of about 1 GeV. 
Unlike the fermionic fluctuations in quark-meson mod¬ 
els, the baryonic fluctuations of the parity-doublet model 
alone, between the UV cutoff and too, are thus not strong 
enough to drive the system into the broken phase. The 
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mesonic fluctuations dominate the flow at T = // = 0. 
In purely mesonic 0(IV)-models, on the other hand, the 
difference between pole and curvature masses is known 
to be a few percent effect and hence negligible [35] [52]. 
We therefore neglect this discrepancy here as well and 
use the standard curvature masses to fix the parameters 
in our calculations. In addition, we have verified, how¬ 
ever, that a larger value for the pion curvature mass, as 
a simple fix to compensate missing wave-function renor¬ 
malization in the vacuum, does not change our results in 
any substantial way. 

From the effective potential in the IR one easily ex¬ 
tracts meson curvature masses as 

m a = ^2U' k + 417"</> 2 , m 7r = y // 2t^, 0. (27) 

The parameters in the UV potential Vm are adjusted to 
realize the physical pion mass = 138 MeV in the IR, 
which essentially determines the parameter c. Due to its 
nature as a broad resonance the mass of the sigma meson 
is not so well constrained, we fix it to reasonable values 
of m a ~ 500 MeV. 

As usual, the Goldberger-Treiman relation is used to 
connect the minimum of the mesonic potential in the 
vacuum at 7 ? = 0, a = a 0 to the pion decay constant, 
do = f-K = 93 MeV. The Yukawa couplings hi and ft .2 are 
then fixed from Eq. ( |16[ ) for a given mo by the masses 
m± of the nucleon (939 MeV) and its parity partner 
(1535 MeV). 

For determining the in-medium condensate EEU5K53] 
at the phase transition one uses the Feynman-Hellmann 
theorem in combination with the Gell-Mann-Oakes- 
Renner relation. This provides a connection between 
the saturation density of nuclear matter no and the in¬ 
medium chiral condensate o'(no) 


crjrio) 




(28) 


Of course, the ratio a ^°' > also depends on the value of the 
so-called nucleon sigma term a jv- The remaining model 
parameters mo and h v can be used to fix the right satura¬ 
tion density of symmetric nuclear matter no ~ 0.16 fm -3 
and a phenomenologically reasonable value for the nu¬ 
cleon sigma term er/v — 36 MeV [53:. 


Due to the fermionic minus sign, the fermionic contri¬ 
bution to the flow generally generates a negative mesonic 
mass term and hence drives the cr-field away from zero 
and into the broken phase. The integrated fermionic 
flow in the parity-doublet model thereby decreases more 
rapidly with a as compared to quark-meson model cal¬ 
culations. This is because the baryon masses start out 
with a rather large chirally invariant mass too for a = 0 
already and further increase with <r, cf. Eq. (16). Starting 
our flow at A = 1 GeV the baryonic fluctuations therefore 
get suppressed rapidly with <j which results in a compar¬ 
atively large negative contribution to the curvature of the 
effective potential near the origin in field space at cr = 0. 
With the infrared minimum fixed at a o = f -n it then turns 
out that we have to include a small Xo4> 6 term in order to 
generate an overall infrared potential that is sufficiently 
shallow around do, for a sufficiently low curvature mass 
of the a meson around 500 MeV. This problem can be 
avoided by starting the flow at a larger UV scale A. We 
have verified that the value of A6 for constant infrared pa¬ 
rameters indeed decreases with increasing A. The same is 
true for the UV potential in the microscopic bare action 
when we include the mesonic fluctuations which coun¬ 
teract the fermionic ones. In the parity-doublet model, 
where the latter are not so strong, we start at A = 1 GeV 
just in the broken phase already, i.e. with a small nega¬ 
tive mass term in the UV potential as mentioned above. 
For a meson masses around 500 MeV we then occasion¬ 
ally also need a small irrelevant coupling A6 in the UV 
potential, if A is not sufficiently large relative to TOq. 



H B in MeV 


III. RESULTS 

A. Extended mean-field 

We first discuss our results from the extended mean- 
field (eMF) approximation which amounts to only taking 
the fermionic contributions to the FRG flow into account. 
This is comparable to the mean-field approximation pro¬ 
vided that the contributions from the vacuum term are 
properly included in the grand potential [36]. 


FIG. 1. (Color online) The chiral condensate for mo = 
820 MeV as a function of baryon chemical potential /xs. The 
plot focusses on the nuclear-matter transition and shows re¬ 
sults for different vector couplings h v . Increasing values of 
h v shift the transition to larger fis- The binding energy per 
nucleon (vertical dotted line), the in-medium condensate (hor¬ 
izontal dotted line) and the nuclear saturation density of sym¬ 
metric nuclear matter are reproduced for h v = 4.1. 

Our eMF results for the zero-temperature chiral con¬ 
densate d over the baryon chemical potential /i# in 
the region of the nuclear-matter transition with toq = 
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820 MeV and different vector couplings h v are shown in 
Fig- 0 The parameter set with h v = 4.10 reproduces 
the nuclear-matter binding energy of Ef, ~ 16 MeV per 
nucleon (by the location of the discontinuity at /i c B = 
(939 — 16) MeV= 923 MeV marked by the dotted red 
vertical line) together with a saturation density of no — 
0.16 fm~ 3 , an in-medium condensate er(no) — 69 MeV 
(marked by the dotted red horizontal line), correspond¬ 
ing to a nucleon-sigma term of ajy — 36 MeV, cf. Eq. (281. 




FIG. 3. (Color online) Baryon number density as a function 
of baryon chemical potential fis for mo = 820 MeV and h v = 
4.10. 


hardly be relevant for the equation of state of nuclear 
matter in neutron stars. Mesonic fluctuations might well 
change this, however, as they tend to bring both tran¬ 
sitions much closer together as we will see in the next 
subsection. 


FIG. 2. (Color online) Solutions Afis = \h v ujo of the cu-gap 
equation at mo = 820 MeV, Hb = 892 MeV as a function of 
the chiral condensate for different values of the vector coupling 

h v • 

The corresponding solutions for A/xb = ih v Cuo of the 
w-gap equation, cf. App. [A] as a function of a are shown 
in Fig. [5J These were all obtained at /xb = 892 MeV 
which corresponds to the nuclear-matter transition for 
the smallest vector coupling h v = 0.30. This illustrates 
how the first-order phase transition at n c B gets shifted to 
larger values of /xb> and weakened at the same time, by an 
increasing vector coupling, due to a non-zero shift A/xb in 
the chemical potential. In order to achieve this desired 
effect it is crucial, however, to have a non-trivial solu¬ 
tion with non-zero A/xb = ih v Cdo for a values larger than 
the in-medium condensate at saturation density, a(n o), 
which is between 50 MeV and 69 MeV for the parame¬ 
ters used here (cf. Fig. [lj. We also note that the bifur¬ 
cation point in the w-gap equation stays put at around 
<7 ~ 80 MeV for all values of the vector coupling which 
implies that fi c B can not be shifted any further, once a (no) 
reaches this value, because A/xb(^(^o)) remains at zero 
regardless of the size of the vector coupling from then on. 

Fig. [3]shows the baryon density as a function of /xb for 
mo = 820 MeV and h v = 4.10, i.e. with the parameters 
for a phenomenologically correct eMF description of the 
nuclear-matter transition. At fi c B = 923 MeV the den¬ 
sity jumps discontinuously from zero to no — 0.16 fm _ , 
the saturation density of symmetric nuclear matter. The 
second discontinuous transition (not shown in Fig. [TJ is 
observed at — 1472 MeV. It occurs at a density of 
around 13.7 times the saturation density and would thus 



FIG. 4. (Color online) The chiral condensate as a function of 
fj,B for the nuclear-matter and for the chiral phase transition 
of the model, with physical pion mass and in the chiral limit. 

That the second transition within the high-density 
phase is indeed the chiral transition is demonstrated in 
Fig. [4] which shows the chiral condensate over the baryon 
chemical potential again: The upper blue curve corre¬ 
sponds to the result of Fig. [l] with correct nuclear-matter 
parameters, the pion mass was thereby adjusted to its 
physical value, m n = 138 MeV. The lower green curve 
was obtained with the same parameters except for the 
explicit chiral symmetry-breaking parameter which was 
set to c = 0 for the chiral limit with ra, = 0. The fact 
that a then drops to zero in the second transition allows 
to unambiguously identify this as the chiral first-order 
phase transition of the model inside the high baryon- 
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H B in MeV H B in MeV 


FIG. 5. (Color online) Masses of the nucleons and their parity 
partners at mo = 820 MeV for physical pion mass and in the 
chiral limit. 


density phase. These results also allow to deduce the 
parameters of the nuclear-matter transition in the chiral 
limit, with <r 0 ~ 80 MeV, Ef, ~ 47 MeV, no — 0.064 fm -3 
and <r(no) — 0.89 <7q. In Fig. [5] the corresponding masses 
of the nucleons and their parity partners as a function 
of hb are shown for both, the physical pion mass and in 
the chiral limit as well. It is evident that they become 
degenerate at large /tb, after the second phase transi¬ 
tion. For a physical pion mass the splitting of the baryon 
masses stays finite after the chiral phase transition and 
then smoothly tends to zero with m± —> mo. This con¬ 
firms the mechanism through which chiral symmetry is 
realized in the parity doublet model. Rather than yield¬ 
ing vanishing masses in the chiral limit the baryon masses 
become degenerate at m± = mo, as observed in the lat¬ 
tice QCD simulations upon chiral symmetry restoration 

[IMS]. 


B. FRG Results with Mesonic Fluctuations 

In this section we present our results obtained from 
the full FRG flow, the solutions of the full flow equation 
discretized on a grid in field space |13| . including the 
fluctuations from collective mesonic excitations as per 
Eq. ((25]). 

In Fig. [6] the resulting expectation value of chiral con¬ 
densate ef is shown as a function of baryon chemical 
potential /r b■ The two representative data sets corre¬ 
spond to two different values for mo = 800 MeV and 
900 MeV with different vector couplings h v = 0.3 and 
h v = 0, respectively. In both sets the other parame¬ 
ters were adjusted to reproduce the physical pion and 
baryon masses. For mo = 800 MeV the nuclear-matter 
and chiral phase transitions occur at — 866 MeV 

and /i c J x ^ — 896 MeV, respectively. For mo = 900 MeV 
these values are ~ 923 MeV for nuclear matter and 


FIG. 6. (Color online) The chiral condensate as a function 
of fiB for the nuclear-matter and the chiral phase transition 
with mo = 800 MeV and parameters fixed to the in-medium 
condensate (horizontal dotted line), and mo = 900 MeV with 
parameters fixed to the binding energy per nucleon (vertical 
dotted line). 

^M,\) ^ 1009 MeV for the chiral transition. The corre¬ 
sponding nuclear-matter saturation densities for the two 
parameter sets are uq ~ 0.10 fm -3 with mo = 800 MeV 
and no — 0.01 fm -3 with mo = 900 MeV. They are both 
much smaller than the phenomenological no — 0.16 fm -3 
which appears to be due to missing density contributions 
from the omega meson. As can be seen in this figure 
we generally observe that one can either reproduce the 
physical binding energy per nucleon (as indicated by the 
dotted red vertical line) or the correct in medium con¬ 
densate (the dotted red horizontal line) but not both at 
the same time. 



g B in MeV 


FIG. 7. (Color online) Masses of the nucleons and their parity 
partners at mo = 800 MeV and 900 MeV with physical pion 
masses, full FRG results with parameters as in Fig. [6] 

The corresponding baryon masses are shown for the 
same two parameter sets in Fig. [7] As before, they both 
tend to mo from above and below beyond the chiral tran- 
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FIG. 8. (Color online) Solutions of the cj-gap equation for different vector couplings h v with mo = 800 MeV (left) and mo = 900 
MeV (right), and values for the chemical potential at the respective zero-temperature nuclear-matter transitions is = 866 MeV 
and fiB = 923 MeV, as in Fig. [6] 


sition which is generally closer to the liquid-gas transition 
than in the eMF results as mentioned above. 

The reason why, in contradistinction to mean-field and 
eMF calculations, the repulsive vector interaction is inef¬ 
fective in shifting the nuclear-matter transition with full 
mesonic fluctuations becomes clear from Figs. [8a|and|8b| 
In these figures we plot the shifts A/r# in the baryon 
chemical potential from the solutions to the w-gap equa¬ 
tions over the chiral condensate for various vector cou¬ 
plings with the two values of mo used in Fig. [6] 

For mo = 900 MeV the bifurcation point in the w-gap 
equation always lies below the value of the in-medium 
condensate cf(no) at the respective saturation density. 
In this case the regime of non-vanishing A= ih v ioo 
does not extend beyond the position of the minimum 
in the cr-direction of the effective potential at satura¬ 
tion density, and larger vector couplings thus have no 
effect on the nuclear-matter transition. In the case with 
mo = 800 MeV, even though the bifurcation point at 
a ~ 75 MeV lies slightly above a(no) — 69 MeV, the 
repulsive vector coupling only has a very minor influ¬ 
ence on the in-medium condensate. For example, with¬ 
out vector coupling we obtain with this parameter set 
a (no) ~ 67 MeV which we can shift to the phenomeno¬ 
logical (J(no) — 69 MeV by increasing the vector coupling 
from h v = 0.0 to h v = 0.3. Larger vector couplings will 
not lead to any substantial further shift in the nuclear- 
matter transition. In addition, we have verified that the 
bifurcation points in the w-gap equation remain where 
they are in Figs. [8a| and |8b| also for much larger values of 
the vector coupling h v than those shown there. We con¬ 
clude that the non-vanishing A/r# values essentially only 
change the form of the local potential away from the min¬ 
imum. Depending on the chiral nucleon-mass parameter 
mo they do either not at all or not significantly influence 
the grand potential which is evaluated at <tq and a(no) 


on either side of the transition. 


C. Finite Temperature Results 


Since the flow equation for the effective potential in 
the parity-doublet model of the previous section is al¬ 
ready formulated for finite temperatures it is in princi¬ 
ple straightforward to obtain the phase diagram of the 
model in the (T, /rs)-plane. In this subsection we present 
first results from the full RG flow at finite temperature 
(and chemical potential). As the impact of the repulsive 
vector-meson interaction can be neglected in full RG cal¬ 
culations, as we have seen, it is not included here, i.e. we 
use h v = 0 in the following. These results are only meant 
to give a rough qualitative estimate of the critical tern- 

(n) (V) 

peratures T c for the nuclear-matter and T c for the 

chiral transition. We have chosen the two parameter sets 
from Sec. Ill B| (cf. Fig. [6]), both with h v = 0 here, to 
perform the finite temperature calculations presented in 
this subsection. The results are summarized in Figs. [9a| 
and |9b| as well as in Figs. |10a| and |10b[ 

In Fig. [9a] the chiral condensate for mo = 800 MeV 
(with h v = 0 and the other parameters as in Sec. Ill B) 
is plotted over the (T, p.^)-plane. As expected, both 
first-order phase transitions, the nuclear matter and the 
chiral transition, get weaker with increasing temper¬ 
ature and eventually become continuous. The corre¬ 
sponding critical lines in the (T, /i b )-plane are shown in 
Fig. 10a The location of the critical endpoint (CEP) of 


the nuclear-matter transition can be estimated to be at 


T, 


(n) 


19 MeV and /i C g ~ 871 MeV. In comparison, 


that of the chiral CEP is around T', 


(x) 


18 MeV and 


^B 


(x) ~ 896 MeV. 


The analogous plots of the chiral condensate and the 
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FIG. 9. (Color online) Chiral condensate over chemical potential fiB and temperature T from the full FRG flow for the two 
different parameter sets with mo = 800 MeV (left) and mo = 900 MeV (right) from Sec. IIIB both with h v = 0. 
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FIG. 10. (Color online) First-order lines with the corresponding critical endpoints in the ( T , p_B)-plane for the liquid-gas (green) 
and chiral (red) phase transitions using the same data sets as in Figs. 9a and 9b for mo = 800 MeV (left) and mo = 900 MeV 
(right). 


critical lines for the parameter set with m 0 = 900 MeV 
are shown in Figs. |9b| and |10b| In this case the nuclear- 
matter transition is much weaker in the first place. It 
thus also turns continuous already at a very low temper¬ 
ature of T c (n) ~ 2 MeV with ~ 925 MeV. We esti¬ 
mate the location of the chiral CEP to now be around 
T c (x) ~ 19 MeV and f/ B {x) ~ 1012 MeV. 

Relative to their corresponding zero temperature tran¬ 
sitions at /ig"' ~ 866 MeV and fi°g X ^ — 895 MeV for 
mo = 800 MeV, or fig"* ~ 923 MeV and — 

1015 MeV for mo = 900 MeV, especially the nuclear- 
matter CEPs tend to appear at somewhat larger hb val- 

(y) 

ues. The chiral first-order transitions from T = 0 to T c 
basically almost follow lines of constant /ig, especially 
for mo = 900 MeV. 


The slope of the first-order lines in the (T, /Lt/j)-plane 
is determined by a Clausius-Clapeyron equation [56j . 


dT c An 
dn c As 


(29) 


For the discontinuity in the number density we have 
An > 0 with increasing fj,B across the first-order lines in 
both cases. For an order-disorder transition one would 
expect the entropy per particle to increase, and with 
An > 0 hence also the discontinuity in the entropy den¬ 
sity to be larger than zero, i.e. As > 0. This is the 
typical behavior of the chiral transition line in mean-field 
studies. Here, As in the chiral transition with mesonic 
fluctuations and explicit symmetry breaking appears to 
be very small. It may well be positive in the chiral limit. 
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More interestingly, because the entropy per particle 
decreases from the gaseous to the liquid phase, it is pos¬ 
sible to have As < 0 in the liquid-gas transition of nu¬ 
clear matter despite the fact that the number density in¬ 
creases across the transition. This is what we obtain for 
the nuclear-matter transition in the parity-doublet model 
with mesonic fluctuations for mo = 900 MeV and at low 
temperatures also for mo = 800 MeV. It is a genuine 
effect of mesonic fluctuations that they tend to change 
the sign in the discontinuity of the entropy density as 
compared to mean-field studies. It is known from the 
quark-meson model where they change the mean-field 
chiral transition into a transition to bound quark mat¬ 
ter m- The same effect turns the relativistic analogue 
of a Chandrasekhar-Clogston transition inside the pion 
condensation phase at finite isospin chemical potential, 
as observed at mean-field with As > 0, into a first-order 
transition to a stable Sarma phase with As < 0, when 
mesonic fluctuations are included EZIES]. This would 
be analogous to a partially polarized phase in ultracold 
Fermi gases at unitarity. For the nuclear-matter tran¬ 
sition in QCD it appears to be rather unusual to have 
As < 0. It is not observed in chiral effective field the¬ 
ory for example [59]. Since An tends to be too small 
in our FRG results with mesonic fluctuations, however, 
especially for mo = 900 MeV, the slope of the nuclear- 
matter transition line might still turn out to be negative 
with fluctuations as well, if a more realistic An, with 
An ~ 0.16 fm -3 at T = 0, is sufficient for the entropy 
density to also increase across the transition (leading to 
As > 0) despite the fact that the entropy per particle 
should decrease. 


IV. CONCLUSIONS 

We conclude, that the inclusion of a heavy parity part¬ 
ner in a chiral baryon-meson model such as the parity- 
doublet model within an FRG framework allows for a si¬ 
multaneous description of the liquid-gas transition of nu¬ 
clear matter together with a chiral first order transition 
inside the high baryon-density phase. The quantitative 
properties of symmetric nuclear matter are well repro¬ 
duced in the extended mean-field approximation without 
collective mesonic fluctuations. Including mesonic fluc¬ 
tuations does not change the qualitative conclusion of 
the existence of two distinct first-order phase transitions 
which can be identified as a liquid-gas transition of nu¬ 
clear matter and the chiral phase transition at which the 
nucleons become degenerate with their parity partners. 

First calculations at finite temperature (and chemical 
potential) provide the general features of the phase dia¬ 
gram of the model with the two first-order lines ending in 
two distinct critical endpoints. The inclusion of mesonic 
fluctuations thereby has effects, especially on the liquid- 
gas transition, that are known from quark-meson models. 
As compared to (extended) mean-field studies they lead 
to a sign-change in the slope of the critical line indicating 


a sign-change in the discontinuity of the entropy density. 

In contrast to the case with purely baryonic fluctu¬ 
ations, however, the repulsive iso-scalar vector meson 
interaction turns out to be inefficient in adjusting the 
binding energy per nucleon of symmetric nuclear mat¬ 
ter. As a result, the FRG treatment of the parity-doublet 
model with mesonic fluctuations does not provide a quan¬ 
titatively fully successful phenomenological description 
of the nuclear-matter transition at this point. A pos¬ 
sible remedy which has proven to work for the chiral 
Walecka model [SI IB] would be to parametrize the ef¬ 
fective potential in the vacuum and to consider only ther¬ 
mal fluctuations. In a longer term one might embed 
the parity-doublet model into a (Polyakov-)quark-meson- 
baryon model to include the effects of fluctuating light 
quark degrees of freedom in the chirally restored and de- 
confined phase and at the initial microscopic scales. 
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Appendix A: Gap Equations 


For the grand potential we need to minimize the re¬ 
sulting effective action with respect to the cr-meson field. 
The w-meson is not included as a fluctuating field here, 
but treated in a stationary phase approximation with a 
complex saddle point for the wo-integration as mentioned 
in the main text. For both, the minimum of the effective 
potential in the cr-direction and the (purely imaginary) 
saddle point in ui o, we thus require the partial derivatives 
of Uq, where Uq = Uk for A: —>• 0, with respect to a and w>o 
to vanish when a = a and wo = Wo- The corresponding 
two gap equations are of the form, 


j- ( / dkd k U k 


a—a,LOQ—ujQ 


and 


d_ 


[ dkdkU k 


= o ■ 


(Al) 


(A2) 


a=cr ,ujq—ujq 


Vice versa Eqs. (A2| and (All define the expectation 
values of the meson fields a and Coq. 

The gap equation for ujq can be derived from the 
fermionic RG flow in Eq. (24). The derivative with re¬ 


spect to wo is thereby equivalent to a derivative with 
respect to the baryon chemical potential 


d 

duj 0 



(A3) 
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For T —> 0 the gap equation for to o can be obtained ana¬ 
lytically in closed form from the expression, 


d 

du } 0 


dkd k U k 


= ih 4 fc^ 4 7 ± (fcf) 

"67T 2 A. | 4 e±(fc)| fe=fe ±| 


with 


±,n _ 2(fc 2 +?n 0 -ej ) + cr 2 (/i? + /il) 
7()_ ^(ef-ef) 


(A4) 

(A5) 


Here, kp are the Fermi momenta defined by 
e fe (^f) = &B- 


(A6) 


Eq. (A4) also allows us to obtain the baryon number 
density, 


n B = — 


dU 0 

d^B 


1 

i h v dujp 


dkd k U k 


(A7) 


With for the contributions to the number densities 
of nucleons and their parity partners separately we thus 
find, 


n B 



4 

6 71-2 ± I4efe(fc)| fe=fc ±| ' 


(A8) 
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